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Abstract 

We consider a lattice of coupled circle maps, a model arising natu- 
rally in descriptions of solid state phenomena such as Josephson junc- 
tion arrays. We find that the onset of spatiotemporal intermittency 
(STI) in this system is analogous to directed percolation (DP), with 
the transition being to an unique absorbing state for low nonlinearities, 
and to weakly chaotic absorbing states for high nonlinearities. We find 
that the complete set of static exponents and spreading exponents at 
all critical points match those of DP very convincingly. Further, hy- 
perscaling relations are fulfilled, leading to independent controls and 
consistency checks of the values of all the critical exponents. These 
results lend strong support to the conjecture that the onset of STI in 
deterministic models belongs to the DP universality class. 
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1 Introduction 



Spatiotemporal intermittency (STI) in coupled map lattices (CML) has been 
extensively studied in varied contexts, especially as it is the precursor of fully 
developed spatiotemporal chaos in extended dynamical systems [1]. There 
are two types of motion seen in systems exhibiting STI: laminar and turbu- 
lent. The laminar region is characterized by periodic or even weakly chaotic 
dynamics, while no spatiotemporally regular structure can be seen in the 
turbulent regime. A laminar or 'inactive' site becomes turbulent or 'active' 
at a particular time only if at least one of its neighbours was turbulent at an 
earlier time, i.e., there is no spontaneous creation of turbulent sites. Hence 
a turbulent site can either relax spontaneously to its laminar state or con- 
taminate its neighbours [2] . This feature is analogous to directed percolation 
(DP). Also, once all the sites relax spontaneously to its laminar state, the 
system gets trapped in this state for all time. Hence, the laminar state is "ab- 
sorbing" in STI. The existence of the absorbing states led to the conjecture by 
Pomeau that STI in deterministic models also belongs to the DP universality 
class [3]. To be more precise, if the propagation rate of turbulence is below 
a certain threshold, the turbulent states die out and the system remains in 
the laminar state for all time (laminar phase/inactive phase). On the other 
hand, on exceeding this threshold, the turbulent states start "percolating" 
in spacetime (turbulent /inactive phase). 

While there is substantial evidence of DP universality in stochastic mod- 
els exhibiting a continuous transition to an absorbing state [4, 5], it is of 
considerable interest to examine the robustness of DP critical behaviour in 
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systems with completely deterministic evolution rules. To this end, Chate 
and Manneville introduced a simple CML [6] exhibiting STI and possessing 
infinitely many absorbing states. Surprisingly, it was found that not only 
were the exponents governing the onset of STI different from those of DP, 
they were also non-universal in nature [6]. This non-universal behaviour was 
considered to be due to the existence of travelling solitary excitations (soli- 
tons) with long life times in this model. It was argued that the existence of 
these sohtons, which spoiled the DP nature of the universality class, was an 
artefact of the synchronous updating rule and asynchronous updates were 
tried with this model. This asynchronous update destroyed the solitons and 
then the static exponents were found to be consistent with DP [7]. Later it 
was also observed that even finite lifetime solitons could completely change 
the nature of transition in a weak soliton region [8] . This indicates that it is 
non-trivial to map deterministic dynamics to stochastic behaviour as various 
spatiotemporal structures (such as these solitons) may introduce long range 
correlations thereby ruining the analogy. Hence it is of considerable interest 
to find CMLs with regimes where there are no additional special spatiotem- 
poral structures, which can be used as clean testbeds for checking the validity 
of the DP universality class. 

Since the state variables of CMLs can often be identified with physi- 
cal quantities (such as voltages, currents, pressures, temperatures, concen- 
trations or velocities) in fairly realistic situations, it is conceivable that 
such models may suggest various experimental possibilities for observing DP, 
which still remains an outstanding problem [4]. Further STI is a common 
phenomena of many extended systems, and is seen for instance in experi- 
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merits on convection [9] and in the "printers instability" [10]. So if the onset 
of STI does in fact exhibit DP universahty, then it could lead to promising 
candidates for observing DP in real phenomena. And indeed a recent experi- 
ment by Rupp et al [11] finds agreement with DP exponents at the transition 
to STI in a 1-dimcnsional system of ferrofluid spikes driven by an external 
oscillating magnetic field. 

Now, in the last decade only the Chate-Manneville CML and its variants 
[8] have been studied in this context. Numerical evidence from more varied 
sources is required, especially in the absence of analytical results, in order to 
settle the question of DP universality in transitions occuring in deterministic 
systems. Thus, it is of considerable interest to be able to find examples of this 
correspondence in systems quite distinct from the Chate-Manneville class, 
and with qualitatively different absorbing states, that would lend credence 
to the Pomeau conjecture. This work provides one such example. 

In this paper, we consider the coupled circle map lattice [12] which has 
been used to model mode-locking behaviour of the type seen in coupled os- 
cillator systems and in diverse experimental systems such as Charge Density 
Waves and Josephson- Junction arrays [13]. We find that this system has 
regimes which show spatio-temporal intermittency, with both unique and 
weakly chaotic laminar regions, at different values of the nonlinearity param- 
eter. Importantly, the solitons which spoilt the DP behaviour in the earlier 
studies are completely absent here. Thus we have at hand a CML without 
any potentially problematic coherent spatiotemporal structures, showing the 
onset of spatiotemporal intermittency very cleanly. This CML can then serve 
as a good testbed for checking the validity of the DP universality class for 
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both unique and fluctuating absorbing states. 

We will now present results from this system strongly indicating that 
the complete set of static exponents characterising the transition to STI are 
completely consistent with DP, both for unique and weakly chaotic absorbing 
states. We will also show that the spreading exponents (dynamic exponents) 
for unique as well as fluctuating absorbing states agree within 3 % of those 
obtained in DP. Further we will demonstrate that the hyperscaling relations 
in case of the static as well as the spreading exponents are also satisfied. Thus 
we will provide two distinct examples of clean DP universality in transitions 
to STI, one of which constitutes the first known example of this correspon- 
dence in a CML with an unique absorbing state, and the other constitutes 
the first example of this correspondence when there exists weakly chaotic ab- 
sorbing state. 

2 Static (Bulk) Exponents 

First, recall the coupled circle map lattice[14]: 

9i,t+i = (l-6)/(^,,t) + |(m_i,,) + m+i,t)) niod 1 (1) 

where t is the discrete time index, and i is the site index: i = 1, . . . L, with 
L being the system size. The parameter e gives the strength of the diffusive 
coupling between site i and its two neighbours. The local on-site map is 
given by 

f{e) = ^ + a; - sin(27r^) (2) 
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where the parameter k gives the nonhnearity. This CML has been studied 
extensively with parallel up-dates and has a rich phase diagram with many 
types of attractors and strong sensitivity to initial conditions[14, 15]. In par- 
ticular, this system also has regimes of spatio-temporal intermittency (STI) 
when evolved parallcly with random initial conditions. Figs. 1 and 2 show 
space-time plots of the spatio-temporal intermittency observed in two differ- 
ent STI regimes. It is clear that no travelling wave soliton-like structures are 
seen in this regime. Hence, it is not necessary to introduce any asynchronic- 
ity here to destroy "sohtonic" behaviour, as this system is naturally free of 
such spatiotemporal excitations in the parameter region studied here. 

We shall now study the onset of spatiotemporal intermittency in this sys- 
tem. Interestingly, as mentioned before, two qualitatively distinct absorbing 
regions can be found in this system. 

(i) When nonhnearity parameter k — 1, there are regions of (e — cu) space 
where the system goes to the synchronised spatiotemporal fixed point 6* = 
^sin~^(^). This constitutes an unique absorbing state (see Fig. 1). We 
closely scrutinise the critical behaviour at 2 critical points in this regime: cu — 
0.064, e = 0.63775, and u; = 0.068, e = 0.73277. These mark the transition 
from a laminar phase to STI. The turbulent sites here are those which are 
different from 6*. 

(ii) When nonhnearity parameter k = 3.1, there are regions of (e — uj) 
space where sites with any value less than 1/2 constitute the absorbing states, 
and sites whose values are greater than 1/2 constitute the turbulent states 
(see Figs. 2 a-b). So now the absorbing states are inifinitely many, as also 
weakly chaotic. We study the critical behaviour at 2 critical points in this 
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regime: uj = 0.18, e = 0.701, and u = 0.19, e = 0.65612. 

As mentioned earlier, we initiate the evolution with random initial condi- 
tions and let the system evolve under parallel updates. The DP universality 
class is characterised by a set of critical exponents which describe the scaling 
behaviour of the quantities of physical interest. The physical quantities of 
interest for such systems are (a) the escape time r, which is the number of 
time steps elapsed before the system reaches its laminar state and (b) the 
order parameter, m(e, L, i), which is the fraction of turbulent sites in the 
lattice at time t. Prom finite-size scaling arguments, it is expected that r 
depends on L such that 

log L laminar phase 
t{u;,€) — <^ critical phase 

exp turbulent phase 

Here, c is a constant of order unity, and the critical point is identified as 
the set of parameter values at which r shows power-law behaviour with z 
being the critical exponent. At Ec the critical value of the parameter e (other 
parameters being held fixed), the order parameter m{e,L,t) scales as 

m ~ (e - CcY, e -> e+. (3) 

when the critical line is approached from above. Also, at the critical ec, the 
order parameter is expected to satisfy the scaling relation 

m - (e - ecf, e e+. (4) 

We compute the above quantities for our CML averaged over an ensemble 
of 10^ initial conditions. The dependence of r on L for different values of 



e is shown on a log-log plot in Fig. 3. Fig. 3 shows this dependence at 
the parameter values k — 1, cu — 0.068 (i.e. at parameter values which 
correspond to a unique absorbing state). It is clear from the graph that an 
algebraic increase can be seen at the value €c — 0.63775. A similar analysis 
was carried out for the parameter values k = 1, uj = 0.064, where a unique 
absorbing state can again be seen and gave the critical value Cc = 0.73277. 
Weakly chaotic absorbing states were seen at the parameter values k — 3.1, 
a; = 0.18, and a; = 0.19. Here the critical values of parameters turned out to 
be €c — 0.70100 for the first case and €c — 0.65612 for the second case. The 
critical exponent z was estimated at these critical values for all four cases. 
The log-log plot of the escape time r against the system size L is shown in 
Fig. 4. It is clear that the same exponent z is seen for all four cases and 
turns out to lie in the range 1.58 — 1.59 which is completely consistent with 
the DP value for this exponent. (See Table I). 

The order parameter is expected to obey the scaling relation 

m(ee, L, t) « r^/'^^ (5) 

for i << T. Therefore the log- log plots of m as a function of time t for 
various lattice sizes must fall on one line when t « t and the power —P/vz 
must correspond to the slope of the graph for these regimes. The order 
parameter is plotted as a function of t on log-log plot in Figs. 5 and 6, where 
the data in Fig. 5 is obtained for the parameter values k — l,uj — 0.068 and 
= 0.63775 (the unique absorbing state case, and that in Fig. 6 is obtained 
for k = 3.1, to = 0.19 and ec = 0.65612 (the case with weakly chaotic 
absorbing states). In both cases, the data for L = 50,100,300,500,1000 
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collapses on to one line, the slope of which gives —P/vz = —0.16. 

The order parameter of systems which belong to the directed percolation 
universality class satisfies the scaling function 

m{e,,L,t)r^L<'/''gUt/L^). (6) 

at the critical value e = €c. We plot the order parameter for our CML at 
the critical values above in Figs. 7 and 8. The data with scaled variables 
M = mL^ and T = t/L^ fall on one curve for various lattice sizes indicating 
dynamical scaling (see Figs. 7 and 8, parameter values are as given in the 
figure captions). This further substantiates the claim that the behaviour of 
the sine circle map CML falls in the directed percolation universahty class. 

The exponent u can be extracted independently, by using the scaling 
relation 

r(L,5)«07(L/0), (7) 

where is the correlation length which diverges as ~ and 5 is given by 
e — ec [16]. Therefore, v can be obtained by adjusting it's value till the scaled 
variables L5'^ and t5^^ collapse onto a single curve. Thus, the exponent (3 
can be obtained from equation 6. 

To extract further critical exponents, we obtain the correlations from the 
pair correlation function given by: 

1 ^ 
^ i=i 

where the brackets denote the averaging over different initial conditions. At 
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criticality one expects an algebraic decay of correlation, i.e. 

where rj' is the associated critical exponent. The log-log plot of the spatial 
correlation function at k = 1, u = 0.064, ec = 0.73277 can be seen in Fig. 9. 
The log-log plot of the correlation function approaches a straight line with 
slope 1 — 77' at large times. The value of the exponent Tyturns out to be 0.302, 
which is consistent with the directed percolation value. The values of the 
exponent at the other critical set of parameter values are hsted in Table I. 
Directed percolation like behaviour is observed for the entire set. 

Thus we have obtained the complete set of static (bulk) exponents, namely 
z, /3, u, rj', that characterize the DP class ( see Table I ). Clearly the values 
of the exponents obtained for the coupled circle map lattice are in excellent 
agreement with the DP values for both unique as well as the weakly chaotic 
absorbing states. 

The exponents also satisfy the hyperscaling relation, 2/?/i/ = d — 2 + rj' , 
where d — 1. 

3 Spreading Exponents 

In the previous section we obtained the static exponents, also known as the 
bulk exponents for our CML, and found that they were in good agreement 
with those of DP. We shall now compute a set of dynamical exponents, called 
the spreading exponents, from the temporal evolution of a nearly absorbing 
system with a localized disturbance i.e. with only few contiguous active 
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0.068 


0.63775 


1.580 


0.28 


1.10 


1.49 


1 


0.064 


0.73277 


1.591 


0.28 


1.10 


1.50 


3.1 


0.18 


0.70100 


1.597 


0.26 


1.12 


1.50 


3.1 


0.19 


0.65612 


1.591 


0.28 


1.10 


1.49 


DP 






1.58 


0.28 


1.10 


1.51 



Table 1: Critical static exponents of the synchronously updated coupled cir- 
cle map lattice for 4 critical points. The first two critical points correspond 
to transitions to an unique absorbing case, while the third and fourth points 
correspond to weakly chaotic absorbing states. The last row shows the cor- 
responding exponents of directed percolation. 
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(turbulent) sites in an otherwise absorbing state. The quantities of interest 
are, the time dependence of N{t), the number of active sites at time t av- 
eraged over all runs, P{t), the survival probability, or the fraction of initial 
conditions which show a non-zero number of active sites (or a propagating 
disturbance) at time t and R'^{t), the mean squared deviation from the origin 
of the turbulent activity averaged over surviving runs alone. The spreading 
exponents are obtained from the time dependence of these quantities which 
show scaling behaviour at criticality. At criticality, we have, 

N{t) f« t", P{t) ^ t-\ B?{t) ^ 

Also, S = P/uz. We shall now compute these quantities for our system and 
compare them with those of DP. 

For A; = 1, we have already seen that the absorbing state is the synchro- 
nised spatiotemporal fixed point. However, when we start with a lattice with 
only a single active site and all other sites at the spatiotemporal fixed point, 
we see that the sytem goes to its absorbing state in about 10 time steps of 
evolution. Further the symmetric diffusive coupling in the CML indicates 
that the temporal spreading will be symmetric about the single active site, 
which is not desirable. Therefore, for good statistics and to counter the sym- 
metric spreading, we need at least two contiguous active sites (see Fig. 10). 
So in our calculations, we have started with two or more contiguous active 
sites, while the background is fixed at 9*. We find that the full set of spread- 
ing exponents obtained thus (see Figs. 12-14) agree within 3 % of the DP 
values (see Table II). 

For A; = 3.1 there is no unique absorbing configuration as above, and 
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any initial configuration with values less than 1/2 constitutes an absorbing 
background. We find that we can now obtain reasonable asymmetric spread- 
ing (and consequently the exponents) with just a single active seed [17], and 
these exponents are the same as those obtained from spreading from two or 
three contiguous active sites (see Fig. 11). Again the complete set of spread- 
ing exponents agree with those of DP within 3 % at all critical points (see 
Table II and Figs. 12-14). 

4 Conclusions 

The evaluation of the complete set of static and spreading exponents at the 
onset of spatiotemporal intermittency in coupled circle map lattices, shows 
that this transition clearly falls in the universality class of directed percola- 
tion. All the critical characteristics of directed percolation, such as hyper- 
scaling, are fulfilled, leading to independent controls and consistency checks 
of the values of all the critical exponents. DP exponents are seen at low 
values of nonlinearity for a unique absorbing state and at high values of non- 
linearity for weakly chaotic absorbing states. It is not necessary to introduce 
asynchronicity or an extra dimension to tune out solitonic behaviour since no 
solitonic behaviour is seen in this model. Thus this model constitutes a very 
clean system where DP exponents are very naturally seen. Model studies 
such as these, showing the clean correspondence between the onset of STI 
and directed percolation, could then lead to promising new candidates for 
observing DP in real phenomena. 
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T] 


5 
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0.068 


0.63775 


0.292 


0.153 


1.243 


1 


0.064 


0.73277 


0.302 


0.158 


1.259 


3.1 


0.18 


0.70100 


0.310 


0.157 


1.272 


3.1 


0.19 


0.65612 


0.308 


0.156 


1.251 


DP 






0.313 


0.159 


1.26 



Table 2: Spreading exponents of the synchronously updated coupled circle 
map lattice for 4 critical points. Two active seeds in an absorbing configu- 
ration is used as initial condition. For the first 2 critical points there exists 
an unique absorbing state, while for the third and fourth points one can 
have many different absorbing states and consequently many different initial 
absorbing backgrounds. However we notice that the exponents obtained are 
quite the same for different initial preparations and thus appears universal for 
2 active seeds. The last row shows the corresponding exponents of directed 
percolation. 
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Figure 1: STI in a synchronously updated coupled circle map lattice of size 
L = 100 with parameters A; = 1, a; = 0.068, €c = 0.73277. The horizontal 

axis is the site index i = 1, . . . L and the vertical axis is discrete time t. The 
absorbing region (black) has sites at the spatiotemporal fixed point of the 
system. 
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Figure 2: STI in a synchronously updated coupled circle map lattice of size 
L — 100 with parameters k — 3.1, u — 0.18, ec = 0.701, where the absorbing 
region has sites below 1/2 and is not unique. The horizontal axis is the site 
index i = 1, . . . L and the vertical axis is discrete time t. The top figure is 
obtained from a density plot of the actual 6 values (the absorbing regions 




Figure 3: Log-log plot (base e) of escape time r versus lattice size L for 
parameters k = 1, to = 0.068 and (a) e = 0.639, (b)e = 0.638, (c) ec — 
0.63775, (d) e = 0.637 and (e) e = 0.635. 




Figure 4: Log-log plot (base e) of escape time r versus lattice size L for all 4 
critical points: (a) = 1, a; = 0.068, = 0.63775 and (b) A; = 1, a; = 0.064, 
Cc = 0.73277 (when there exists an unique absorbing state); (c) k = 3.1, 
u = 0.18, Sc = 0.70100 and (d) k = 3.1, uj = 0.19, = 0.65612 (with 
fluctuating absorbing states) . Table I gives the exponent z of the power law 
fits for the different critical points. 
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Figure 5: Log- log plot of order parameter ■m{t) versus t for different lattice 
sizes at the critical point: k — 1, u; — 0.068 and €c — 0.63775 (when there 
exists an unique absorbing state). For t « r, the data collapses for L ~ 
50, 100, 300, 500, 1000 on to one line, the slope of which gives —l3/yz = —0.16 
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Figure 6: Log- log plot of order parameter ■m{t) versus t for different lattice 
sizes at the critical point: A; = 3.1, a; = 0.19 and €c — 0.65612 (a case with 
inifinitely many absorbing states). For t « r, the data collapses for L — 
50, 100, 300, 500, 1000 on to one line, the slope of which gives —P/vz = —0.16. 
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Figure 7: Log- log plot oi M — mL^/^ vs T — t/L^ a,t the critical point: 
k — 1, cu — 0.068 and Cc — 0.63775 (when there exists an unique absorbing 
state). This rescahng of the order parameter according to Eqn. 6 yields 
independent estimates for P/u and z. The data for system sizes L ranging 
from 2'^ to 2^° collapses onto one curve. 
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Figure 8: Log- log plot oi M — mL^/^ vs T — t/L^ at the critical point: 
k — 3.1, (J — 0.19 and €c — 0.65612 (when there are weakly chaotic absorbing 
states). This rescahng of the order parameter according to Eqn. 6 yields 
independent estimates for P/u and z. The data for system sizes L ranging 
from 2'^ to 2^° collapses onto one curve. 
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Figure 9: Log-log plot (base e) of the spatial correlation function Cj{t) vs j at 
various times t 100, 300, 500, 1000) at the critical point: k = 1, u = 0.064, 
Cc = 0.73277. Cj(t) approaches a straight line with slope 1 — rj' for large 
times indicating an algebraic decay. 
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Figure 10: Spreading of 2 active seeds in an otherwise absorbing lattice of 
size L — 100 at critical point k — 1, u; — 0.068, ec = 0.73277. The horizontal 
axis is the site index i — 1, . . .L and the vertical axis is discrete time t. The 
absorbing region (black) is unique here, with all sites at the spatiotemporal 
fixed point of the system. 




Figure 11: Spreading of 2 active seeds in an otherwise absorbing lattice of 
size L = 100 at critical point k = 3.1, u = 0.18, ec = 0.701. The horizontal 
axis is the site index i = 1, . . . L and the vertical axis is discrete time t. The 
absorbing region (dark) has all sites below 1/2 and is not unique. 
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Figure 12: Log-log plot (base e) of N{t) vs t for all 4 critical points: (a) 
k ^ 1, u ^ 0.068, €c = 0.63775; (h) k = 1, u = 0.064, = 0.73277; (c) 
k = 3.1, a; = 0.18, = 0.70100 and (d) k = 3.1, w = 0.19, = 0.65612. 
Table II gives the exponent r} of the power law fits for the different critical 
points. 
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Figure 13: Log-log plot (base e) of of P{t) vs t for all 4 critical points: (a) 
A; = 1, a; = 0.068, = 0.63775; {h) k = I, u = 0.064, = 0.73277; (c) 
k = 3.1, cu = 0.18, ec = 0.70100 and (d) k = 3.1, u = 0.19, = 0.65612. 
Table II gives the exponent S of the power law fits for the different critical 
points. 
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Figure 14: Log- log plot (base e) of R^{t) vs t for all 4 critical points: (a) 
k ^ 1, u ^ 0.068, ec = 0.63775; (h) k = 1, u = 0.064, = 0.73277; (c) 
k = 3.1, cu = 0.18, ec = 0.70100 and (d) k = 3.1, u = 0.19, = 0.65612. 
Table II gives the exponent Zs of the power law fits for the different critical 
points. 



